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ABSTRACT: We present a unified algebraic Bethe ansatz for open vertex models 
which are associated with the non-exceptional A^n, ^2n-it Lie algebras. By 

the method, we solve these models with the trivial K matrix and find that our results agree 
with that obtained by analytical Bethe ansatz. We also solve the , , D^> models 
with some non-trivial diagonal K-matrices (one free parameter case) by the algebraic 
Bethe ansatz. 
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1 Introduction 

The algebraic Bethe ansatz (ABA) [1]- [3] has been proven to be a more powerful mathe- 
matical method in constructing and solving integrable models. For multi-particles models, 
other methods such as the co-ordinate Bethe anzatz[4], become too cumbersome, partic- 
ulary for open boundary conditions. One important feature of the ABA method is that 
it permits us to present an algebraic formulation of the Bethe states. The construction 
of exact eigenvectors, besides being an interesting problem on its own, is certainly an 
important step in the program of solving integrable systems. This step, however, depends 
much on our ability to disentangle the Yang-Baxter algebra in terms of appropriate com- 
mutation relations. The simplest structure of commutation rules has been discovered in 
the context of six- vertex model[l],[3] and its multistate generalization [5], [6]. 

However, there are many models such as B n ,C n , D n vertex models whose eigenvectors 
are rather complicated to be constructed. Fortunately, Martins and Ramos successfully 
generalized the ABA and threw light on how to construct such general muti-particle 
eigenvectors. They presented a unified algebraic Bethe ansatz for a large family of vertex 
models with periodic boundary and solved them in the way of the ABA [7]- [11]. 

In the framework of Martins's work, recently, the method has been applied to some 
nineteen-vertex like models and Hubbard-like models models with open boundary condi- 
tions. In refs.[12]-[17], the eigenvectors are explicitly constructed and the algebraic Bethe 
ansatz is demonstrated. The results show that the ABA could suit some systems with 
higher rank algebra symmetry. 

There are many models already solved by the so called analytical Bethe ansatz which 
was firstly proposed by Reshetikhin for close chains[18], and then generalized by Nepomechie 
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and his collaborators into the quantum-algebra-invariant open chains [19]. Later, a set 

(2) (2) 

of open vertex models associated with non-exceptional Lie algebras such as A 2r [, A^-i, 
C£\ vertex models[20] were solved by the analytical Bethe ansatz. Compared 
with the analytical Bethe ansatz, the ABA can provide us more rigorous results. So it is 
well worth reconsidering these models by the ABA. 

In this paper we will consider the A^, A^-i, B^\ C£\ D$ vertex models with the 
trivial K matrix and C£\ models with some non-trivial diagonal K-matrices(one 
free parameter case). These models with period boundary conditions have been solved 
by Reshetikhin in terms of analytical Bethe ansatz and by Lima-santos with the help of 
the ABA[21], respectively. Their quantum-algebra-invariant open chains (corresponding 
to the trivial diagonal reflecting K matrix) have already been considered by Nepomechie 
et al.[20]. Here, we will present a unified ABA for these models. We find that the funda- 
mental commutation relations rules have a common form in terms of the corresponding 
Boltzmann weights. As a consequence, the derivations of the eigenvectors, the eigenvalues 
and the associated Bethe ansatz equations also have a quite general character for these 
vertex models. Our results agree with that obtained by the analytical Bethe ansatz. 

We organize our paper as following. In section 2 we demonstrate the ABA for the 
open A^li A^-i, B^\ vertex models and present the solutions to these models 

with the trivial K matrix. The results for B^ , , models with one free param- 
eter diagonal K-matrices are also obtained. A brief summary and discussion about our 
results are included in section 3. Some coefficients and detail derivations are given as the 
Appendices. 

2 The algebraic Bethe ansatz 
2.1 The R matrices and K matrices 

The R matrices for the A^,A^_ X , Bg\C£\DW models [22], [23] can be expressed as [20] 
(u) = a n (u) Y E u ® En + b n (u) ]T E u <g> E j3 + ( Y c n (u, + E ^(u, ij)E a <g> E- u 

+ ( Y d n (u,i,j)+ Y d n (u, i, jfjEij <g> E-q + e n (u) Y Eu <g> Eji 

+6 i if n (u)E ii ®E Ti +(g n (u) Y +g n {u) Y )E ij ®E ji . (1) 

i<j,j¥=i i>3,j¥=i 

The summations run over 1 to q, i + i — q + 1, q — 2n + 1 for A^, B^ and q = 2n for 
A^-i, (Eij) kl = Sit-Sji. The R matrices satisfy the following properties 

regularity : Rn{0) = p n {0)^Vi 2 , 

unitarity : {u)R^ {—u) = p n (u), 

PT- symmetry : V 12 R^(u)V 12 = [R$] tlta (u), 

crossing - unitarity : M 1 (n)R&\u)*Mi(n)~ 1 R$(-u-2ti n ) tl =p n (u + ti n ). 
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Where p n (u) = a n (u)a n (-u), V = <g> E kh Vf x = 8 u 5 jk , $ n = -y^ln - 2kt) for 

42,42-1 and n = -Ikt) for B^\C£\d£\ k = 2n + l,2ra,2ra - 1, 2w + 2, 2w - 2 
for 42, 42-1, -^n 1 ^, -D^, respectively. £j denotes the transposition in z-th space, 
Mi(ra) = M(n) <g) 1 and 



c%e 4 ( n+1 -^ for42, 

X. 4(n+l/2-i)r, fo^ 2 ) /nr(l) n(l) 



with 



z + i, l<z<2±i 



* = ^ (3) 

h-i ^<^< 9 , for 42, i^),/^) 



2 + i n + 1 < i < 2n. for 42-i, 



(4) 



Here one should notice that 42- i R- matrix we used is U q (C n ) invariant, instead of U q (D n ) 
invariant in refs. [23] , [24] , [25] . These R matrices satisfy the Yang-Baxter equation(YBE) [26] 

itfiV - v)R^(u)rS(v) = R%(v)R$(u)R$(u - v), (5) 

R ( u\u) = RW(u) <g> 1, R$(u) = 1 <g> R^(u) etc., 4? = V 12 RnV 12 - For & N x N square 
lattice, if we can find the K matrices satisfying the so called reflection equations 

R ( u\u - v) K- (u)R%\u + v) K- (v) =K- (v)R$(u + v) K- (u)R&\u - v),(6) 

R$(-u + v) K+ (u) M- 1 (n)B&\-u -v- 2£ n ) M (n) K+ (v) 

=K+ (v) M {n)Rf 2 \-u -v- 2£ n ) M~ l (n) K+ {u)R^{-u + v), (7) 

1 2 

where K± (u) = K±(u) <S> 1, K± (u) — 1 <S> K±(u), then the transfer matrix defined as 

t(u) =trK + (u)U(u) (8) 

constitutes an one-parameter commutative family, i.e. [t(u),t(v)] = 0. Here 

U(u) = T(u)K_(u)T-\-u), (9) 
T(u) = R&\u)R$...Rti>(u). (10) 

As indicated by Sklynanin [3] , the transfer matrix is related to the Hamiltonian of quantum 
chain with the nearest-neighbour interaction and boundary terms through the following 
relation 

t'(0) = c c [2HtrK + (0) + trK' + (0)] (11) 
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and the corresponding integrable open chain Hamiltonian takes the form 



ti 2c f 1 trtf + (0) 



(12) 



with Hk t k + i = Vk,k+iR'kk+i( u )\u=o- c c 1S a constant defined by i^-(0) = • zd. 

From Eq.(6) and Eq.(7), we can see that, given a solution K-{u) of Eq.(6), the matrix 



K+(u) = K_(-u-'d) t M(n) 



(13) 



satisfies Eq.(7). The general solutions to Eq(6) for A { £>, A^_ lt B£\ D$ models have 
been obtained in refs. [24], [25]. For Eq.(6), there are three kinds of diagonal K matrices: 
the trivial K matrix, K matrices without free parameter and K matrices with one free 
parameter. We denote these K matrices as K^\n), K^\u,n,p-) and K^\u, n, £_), 
respectively, where p_ are integer number, are free parameters. K^l\n) is a q x q 
identity matrix and can be written as 



K ( l\n) = I. 



(14) 



For the model with the non-trivial diagonal K-matrices have been solved in ref.[16] 
and the A2I-1 R matrix we used is different from that in refs. [24], [25], therefore, we do 
not consider this two models with non-trivial diagonal K-matrices in this paper. We only 
present the non-trivial diagonal K-matrices obtained in refs. [24], [25] for the B£\C£\ D^> 
models as below, 



K {2 \u,n,p^) t = { 



e- u smh({2p_-n-n c )r]-^), (1 < i < p_) 
sinh((2p_ —n — n c )r] + f ), (p~ + 1 < % < q - P-) 
e u sinh((2p_ — n — n c )rj — |). (q + 1 — p_ < i < q) 



(15) 



e n cosh((2p_ - n - n c )rj - |), (1 < % < pJ) 
cosh((2p_ -n-n c )r]+ f), (p_ + 1 < i < q - p_) (16) 
e u cosh((2p_ — n — n c )r) — | ). (q + 1 — p- < i < q) 



where n c — \ for and n c = 2 for C^\D^\ p_ ranges in [2,n] and [2,n — 1] for 
B^\C(p and respectively. For C£\D£\ Eq.(15) do not hold when 2p„ = n + 2. 
The diagonal K matrices with one free parameter for B^ and are 



f e- M sinh(C_ + f)sinh(C--(2n-3)7] + f), (i = 1) 



Al 3) ( U) n,C-)i = <{ 



sinh(C_ - f ) sinh(C- - (2n - 3)r/ + f ) (2 < i < 2n) (17) 
_ e u sinh(C- - f ) sinh(C- - (2n - 3)77 - f ) (i = 2n + 1) 
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and 

e-t sinh(C_ + f), (l<i<n) 
et sinh(C_ - f), (n + 1 < i < 2ra) 
respectively. For -D^ 1 -*, we can get 

' e~ u sinh(C_ + f ) sinh(C- - (2n - 4)r? + f ), (2 = 1) 
< sinh(C_ - |) sinh(C_ - (2ra - 4)r? + f ) (2 < i < 2n) 
k e u sinh(C- - f ) sinh(C_ - (2n - 4)r/ - §), (i = 2ra + 1) 



tfi 3) (u,rc,C-)< 



tfi 3) (u,n,C-)< 



xi 3) ( M ,n,C- 



K?\u,n,(-)i 



e-f sinh(C_ + |), (l<i<n) 

ef sinh(C_ - f), (n+l<i<2n) 

e"f sinh(C_ + f ), (i = 1, 2, • • • , n - 1, n + 1) 

ef sinh(C- - f ) , (i = n, n + 2, • • • , 2n) 

{ e-f sinh(C- + |), (2 = 1) 

ef sinh(C_ - f), (i 

e^sinh(C- + f), (i 
k e^sinh(C_ - f), (i 



' e-f sinh(C- + |), 
ef sinh(C_ - f), 
e^r sinh(C_ + f ), 
k e^sinh(C_ - f), 



2, •••,71- l,n) 
(i = w + l,w + 2, •••,2w- 1) 

(i = 2n) 

(i = l) 
(i = 2, l,n+ 1) 

(z = n, n + 2, • • • , 2n — 1) 

(i = 2n) 



{ e~ u sinh(C- + f ) sinh(C_ + (2n - 4)77 - f ), (i = 1,2,- 
sinh(C_ + f ) sinh(C_ + (2n - 4)r/ + f ), (i = n + 1) 
sinh(C- - f ) sinh(C_ + (2n - 4)t? - f ), (z = 77) 
[ e u sinh(C_ + f ) sinh(C_ + (2n - 4)t? - f ) , (i = n + 2, 



' e~ u sinh(C_ + f ) sinh(C_ + (2t7 - 4)77 - §), (i = 1, 2, • • 
sinh(C_ + f ) sinh(C- + (2n - 4)77 + § ) , (i = n) 
sinh(C_ - f)sinh(C- + (2n -4)t/- f), (7 = 77 + 1) 

, e" sinh(C_ + f ) sinh(C_ + (2n - 4)77 - §) . (i = n + 2, • • 
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Correspondingly, from Eq.(13) and Eqs. (14-25), we can obtain three kinds of diagonal 
K + matrices: K+\n), K+\u,n,p+) and K+\u,n,£+) f° r Eq.(7). Where p + are integer 
number, ( + are free parameters. 



2.2 The vacuum state and commutation relations 

Firstly, we write the double-monodromy matrix (9) as 



U(u) 



( A(u) fli(u) B 2 (u) 

Di(u) A u (u) A 12 (u) 

D 2 {u) A 21 (u) A 22 (u) 

D q _ 2 (u) A q _ 2 i(u) A q _ 22 (u) 

V G{u) d{u) C 2 (u) 



B 2n -! F(u) \ 
A lq _ 2 (u) E^u) 
A 2q _ 2 (u) E 2 {u) 



A 



9-29-2 



'u) E q ^ 2 (u) 



(26) 



•• C,_ 2 (u) A 2 (u) J 
With the help of Eqs. (5,6), we can prove that U(u) in Eq.(26) satisfy the following equation 

R&\u -v)U (m)4i } (« + v) U (v) =U (v)B&\u + v) U {u)R^{u - v). (27) 
Now we Introduce the vacuum state, 



io) = n(i,o,---,o)*, 



(28) 



where (1, 0, • • • , 0) is a 1 x q matrix, t denotes the transposition. 

Applying the double- row monodromy matrix eq.(26) on the vacuum state eq.(28), we 
can find 

£> o (u)|0) = 0, C a («)|0) = 0, G(m)|0> = 0, 
B a (u)\0)^0, E a (u)\0)^0, F(«)|0)^0, 
A aa (u)\0)^0, A ab (u)\0) = (a^b), 

A(u)\0)^0, A 2 (u)|0>^0. (a=l,2,...,g-2) (29) 

From Eq.(29), we can see that D a , C a and B a , E a , F play the role of annihilation operators 
and creation operators on the vacuum state, respectively. The A, A aa , A 2 are diagonal 
operators on the vacuum state. Considering the definition of U(u) Eq.(9), we have 



A(u)\0) = T(u) 11 K_(u) 1 T-\-u) 11 \0), (30) 

+T(u) a+la+1 K4u) a+1 T- 1 (-u) a+la+1 \0), (31) 

A 2 (u)\0) = ^Tiu^K^T-'i-uUO). (32) 
i=i 
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In above equations, the first term of eq.(31) and the previous q — 1 terms of eq.(32) 
can not be calculated directly but it can be worked out by using the following method. 
Taking v = —u in the Yang-Baxter equation, we can get 



T 2 - 1 (- M )i? 12 (2«)T 1 ( W ) = T 1 (u)R 12 (2u)T 2 - 1 (-u). 



(33) 



Taking special indices in this relation and applying both sides of this relation to the 
vacuum state, we find: 



T(u) a+n T 1 (-u) la+1 \0) = f^u) (T 1 (-u)i 1 T{u)i 1 -T{u) a+ i a+1 T 1 (-u) a+la+1 ) |0), 
T{u) q {T-\-u) iq \Q) = M^Uu) (T-'i-u^Tiuh - T{u) qq T-\-u) qq ) |0>, (i ^ 1, q) 



T(u) ql T-\-u) lq \0) = \ fMT-'i-u^Tiu)^ - M^Uu^M^TiuhT-'i-u) 

i=2 



q-1 



+(/i(«)/2(«) E M u ] ~ h(u))T~\-u) qq T(u) 



qq 



i=2 



where 



Ah 



a n {2u) 



h{u) 



c n (2u,q) 



On (2«) 

c n (2u, q)g n {2u) - g n (2u)a n (2u) 



g n (2u)g n (2u) E?=2 M<f } - a n (2u) £g M^R^(2u) 



with arbitrary j G [2, q — 1]. Simple calculations show that 

e u-4 V sin h(2r/) 



h{u) = < 



. -, for 45, C« 

sinn(u — 2(k — l)??) 

e " sinh ( 2 ^) for 4(2) B (i) D {i) 



Introducing two new operators 

A 2 (u) = A 2 {u)-Uu)A{u)-Uu)Y.M^ l) A aa {u), 

a=l 

here we should notice that Mj^~ 1 " > = M^ la+1 , then we have 
A(u)\0) = K4u) 1 [a n (u)] 2N p n (u)- N \0) = wi(u)|0), 



5-2 



(34) 



(35) 



(36) 



(37) 
(38) 



(39) 
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A aa (u)\0) = {K_{u) a+1 - h{u)K.{u){)\b n {u)f N Pn (u)- N \0) = fc-(«) o w( U )|0),(40) 
A 2 (u)\0) = {K_(u) q - Uu) q j2M^(K4u) a+1 - h^K^u),) 

a=l 

-/ 3 H^-Hi}[enH] 2Ar PnH^|0) = w 9 («)|0>. (41) 
In terms of new operators, the transfer matrix (8) can be rewritten as 

t(u) = Wl (u)A(u) + w(u)k+(u)A aa (u) + w q (u)A 2 (u) (42) 

a=l 

with 

tui(u) = + f 3 (u)K + (u) g + f\(u) K + (u) a+1 , 

a=l 

«;(«)#(«) = tf+(«Wi + Mir^H^),, «;,(«) = (43) 

Where the nested k T (u) matrices take the forms K^\n — 1), K^\u,n — l,p T — 1) or 
(u,n — l,( T ) depending on the choice of boundary, u = u — 2rj, ( T = Ct ^ V 01 Ct + V- 
For examples, k^(u) take '{u, n — — 1) and (u, n — 1, ( T ) for Eqs. (15,16) and 
Eqs.(18,20-25), respectively, while it take K^\n - 1) for Eqs.(17,19). 

In order to construct the general m-particle state, we need to find the commutation 
relations between the creation, diagonal and annihilation fields. Rewriting the equation 
(27) in component form 

R(u-)T^u(u)lR(u + )ztu(v)ii 

= U(v)%R(u + )%%U(u)%R(u-)££, (44) 

where the repeated indices sum over 1 to q, u_ = u — v, u + = u + v . In the following 
formulaes, the repeated indices will mean summing over 1 to q — 2 unless there are some 
special notifications. Taking some components of Eq.(44), we can obtain the following 
fundamental commutation relations for those models, 

B a (u)B b (v) + Sabgi(u, v, a)F(u)A(v) + g 2 (u, v, a)F(u)A ab (v) 

= r(u_)£[B d (v)B c (u) + S dc9l (v, «, d)F(v)A(u) + g 2 (v, u, (45) 

A(u)B a (v) = a\(u, v)B a (v)A(u) + a\{u, v)B a (u)A{v) + a\(u, v)B d {u)A da {v) 

+a\(u, v, a)F(u)D a (v) + a\(u, v)F(u)C a (v) + a\(u, v, a)F(v)Da(u)(AQ) 

A ab (u)B c (v) = f( % ):/( !1 .)jB e («)4( !1 ) + ijf( V )ft( ! i)4) 
+R£(u, v) a J b B f {u)A dc {v) + 5 bc Ri(u, v, b)E a {u)A{v) 
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+Ri(u,v,b)E a (u)A- bc (v)+R£(u,v) a jF(u)D f -(v) 
+5 ab R£{u, v)F{u)C c {v) + Rf(u, v) a jF(v)D f -(u) 

+R£(u,v) a cb f F(v)C f (u), (47) 

A 2 (u)B a (v) = a\{u, v)B a (v)A 2 (u) + a 3 2 (u, v)B a (u)A(v) + a 3 {u, v)B d (u)A da (v) 

+al(u, v, a)Ea{u)A(v) + a s 5 (u, v, d)E d (u)A da (v) + a\(u, v, a)F(u)D a (v) 
+a 3 (u, v)F(u)C a (v) + a 3 8 (u, v, a)F(v)D- a (u) + a 3 (u, v)F(v)C a (u), (48) 

A(u)F(v) = b{(u,v)F(v)A(u) + bl(u,v)F(u)A(v) + bl(u,v,d)F(u)A dd (v) 

+b\(u, v)F(u)A 2 (v) + bl(u, v, d)B d (u)B d (v) + b\{u, v)B d (u)E d (v), (49) 

A ab (u)F(v) = bl(u, v)F(v)A ab (u) + 5 ab bj(u, v)F(u)A(v) + R((u, v) d b c a F(u)A dc (v) 

+S ab bl(u, v)F(u)A 2 (v) + R% (u, v)tB d (u)B c {v) + R[(u, v) a d c b B c (u)E d (v) 
+bl(u, v)E a (u)B b (v) + b 2 5 (u, v, b)E a (u)E b (v), (50) 

A 2 (u)F(v) = b 3 (u, v)F(v)A 2 (u) + b 3 (u, v)F{u)A{v) + b 3 (u, v, d)F(u)A dd (v) 

+b 3 {u, v)F{u)A 2 {v) + b 3 (u, v, d)B d (u)B d (v) + b 3 (u, v)B d (u)E d (v) 
+b 3 (u, v)E d {u)B d {v) + b 3 (u, v, d)E d (u)E d (v). (51) 

D a (u)B b (v) = R?(u, v)2B c (v)D d (u) + c\{u, v, a)B- a {v)C b {u) + 5 ab c\{u, v)F{v)G{u) 
+cl(u, v, d)Ba(u)C b {v) + c\{u, v)E a (u)C b {v) + 5 ab c\{u, v)A{v)A{u) 
+c\{u, v)A(v)A ab (u) + 5 ab c){u, v)A(u)A(v) + c\{u, v)A{u)A ab {v) 
+cl(u, v)A ab (u)A(v) + c\ (u, v)A ad (u)A db (v), (52) 

C a {u)B b {v) = RC(u, v)tB d {v)C c {u) + R 2 (u, v)%B c (v)D d (u) + 8 ab c\{u, v, a)F(v)G(u) 
+cj(u, v)B a (u)C b {v) + cj(u, v, a)Ea{u)C b {v) + 5 ab c\(u, v, a)A(v)A(u) 
+R$(u, v) d b c a A(v)A dc (u) + 5 ab cl(u, v, a)A{v)A 2 {u) 
+$ab c K u i v i a)A(u)A(v) + c?(u, v, a)A(u)Aa b (v) 
+R$(u, v) dc a A dc (u)A(v) + R°(u, v)f a A dc (u)A eb (v) 

+S ab cl(u, v, a)A 2 (u)A(v) + cl(u, v, d)A 2 (u)A ab (v), (53) 

B a (u)E b (v) = R\ e {u,v)ZE c {v)B d {u) + R b 2 e (u,v)£B d (v)B c (u) + 5 ab e\(u,v,a)F(v)A(u) 
+R b 3 e (u, v)f a F{v)A dc {u) + 5 ab e\{u, v, a)F(u)A(v) + R\ e (u, v) c b a d F{u)A cd {v) 
+5 ab e\(u,v,a)F(u)A 2 (v), (54) 

where the coefficients Rf(u, v)° d are zeroes except for a = b = c = d,a = c^b = d,a = 
d 7^ b = c and a + b — c + d — q — 1, (a = A, D, etc., i = 1,2, etc.), 

d n (u-,l,d + l)b n (2v) d n (u+,l,a + 1) _ 

g 1 {u,v,a) = — - — — , g 2 {u,v,a) = — — a + a = q-l. 

e n (u„)a n (2v) b n {u + ) 
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The r(u),f(u) and r{u) are given by 

r{u)= ) ^W ^Hu), f{u) = ^ f -R^{u-Ar ] ), f(u) = ^-R^\u), 
e n [u) a n [u) a n [u) e n [u) 

respectively. The other coefficients are omitted here for their long and tedious expressions. 
2.3 The m-particle state 

Inferred from the commutation relation Eq.(45), we can construct the general m-particle 
state for these models as follow. Let 

$^" 6m K • • • , v m ) = B bl ( Vl )^ b r(v 2 , ■■■,v m ) 

m 

+F(vi) E ®t-2 dm (v2, v m )Sf.t^ *}) 

i=2 

xA™" 2 (^; {v u Vi}) gi {v u Vi , b^Afa)^ 

m 

+F(Vl) E $t-2 dm (v2, ■ ■ ■ , Vi, ■ ■ ■ , V m )\f m -\ Vl] {V U ^})t'.c:k lC2 
i=2 

x ^ c n:(^; «i})ff2(ui, v h h), (55) 



where 



j=i+l 



r 



t±_i (« - ^)^:: 1 1 /:: 2 1 (« - ^-0 • • • W (« - «i) (56) 
with s£:::fc( Vl ; v 2 ,---, Vm ) = ut= i <w [^V)U = ai>; «i, ^ • • • , v m ) = n™ i a [(u, Vi ), 

(l = 1,3), $ = 1, ® b i(vi) = B bl (vi). The means missing of v { in the sequence. 
Then the general m-particle state is defined by 

|T m fa, • • • , t/ m )> = KV bm (vu v m )F b -- b -\0), (57) 

which satisfy the n — 1 exchange conditions 

. . . ? Vij . . . ; Vm )F b ^ b - |0> = 
^•••a ia , +1 ... 6m(vi) . . . . . ^^m+i^ _ v . +l)F H-b m \ 0) (58) 

We can prove Eq.(58) by mathematical induction method. 
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2.4 The eigenvalue and Bethe equations 

We can apply the operators x (x = A, A aa , A 2 ) on the eigenstate ansatz and obtain (see 
Appendix B) 

x(«)|T m (?;i, • • • , v m )) = \V x (u, {v m })) 

m 

(u,vf, {v m }) dd ) 

1=1 

m 

/ 2 ) 



(u,vf, {v m }) dd ) 
1=1 

m 

{u,Vi\ {v m }) 

i=l 
m 

+ E U i> a«)l*m-l(M> V i\ { v m})a x aJ 

i=l 

m— 1 m 

i=l j'=i+l 
m— 1 m 

i=l j=i+l 
n—1 m 

U,Vi,Vj] {v m }) dl ) 

i=l j=i+l 
m—1 m 

+ E E ^4,di( u > u <. u j)l*m-2( w > u <> u j;{ u m})di), (59) 
i=l j=i+l 

where the expression of |\^)'s and coefficients -ff^i = 1)2,3,4) are given in Appendix 
B. Using Eq.(42), we then get 

t(u)\T m (v u ■ ■ ■ ,v m )) = ^i(m)cji(m)A™(w; v 1: ■ ■ ■ , v m )\T m (vi, ■ ■ ■ ,v m )) 

+w(u)u(u)A™(u; v u ■ ■ ■ , v m )^- d -( Vl , • • • , v m )n(u; {vm})^^"^) 
+w q (u)u q (u)Af(u] Vi, ■ ■ ■ , v m )\T m (vi, V m )) + wi., (60) 

where u.t. denotes the unwanted terms, 
ri(u-{v m })t: d c : = 

k+(u) a L(u, V2)tt ■ ■ ■ L ^ v m ) h h 2 g l d -k-(u) hm 

xi-^-fi.^^^i-^-u^Oj-^i • ■ -i^C-u,^)^. (61) 
with Vi = Vi — 2i], and 

= *-"[?- (62) 

p n -l{U - V) 
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Thus, we get the conclusion that|T m (wi, • • • ,v m )) is the eigenstate of t(u), i.e. 

t(u)\T m (vi,"-,v m )) = {w 1 (u)u 1 (u)A^(u;v 1 , ■ ■ ■ ,v rn ) 
+w{u)uj(u)A™(u; ui,---, u m )ri(u; {£ m }; {v£J}) 
+w (? (u)u; (? (w)A^(u; v ± , ■ ■ ■ , v m )}\T m (v!, v m )) 
= r(u;{v m })\r m (v 1: ---,v m )), (63) 

if the parameters satisfy 

ri(fi; {v m })F b -- b ™ = T^u; {v m }; {v^})F b ^ b - , (64) 

ri(^;{v m };{<^}) = -p n ^i(0)— — — T7 — 7— jy^i(^)- (« = l,---,"i) (65) 
All unwanted terms cancel out by the following three kinds of identities 

p L ( v ) _ T (di)( v } w i( u ) a l( u > V A + Si! + up ^ 

* wi(u)al(u, Vi) + EL? Wd+i(«)i^(u, t>i)S + w <?( w ) a !K v ' 

^ / V N = r (di) w d 1+ i(u)Rj(u, v u di) + ^(m)o|(m, Uj, di) 
1 Wj 1+1 (u)R£ (u, Vi, di) + w g (u)al(u,Vi, d x ) ' 

q q 
i=i i=i 

-E w iHH^ dl (u, v h Vj)\^{vj) + E^(ttX(«,«i,*(«i)A(fi) = (68) 
z=i ;=i 

with #1 = A, = A u , x q = ^2, ^d+i(' u ) = w(u) (u) , di — 1, 2, • • • , q — 2 do not sum 
in Eqs.(66,67,68). 

From Eqs.(61), (63) and (64) , we can see that the diagonalization of r{u) is reduced 
to finding the eigenvalue of Ti(u;{v m }) which is just the transfer matrix of the same 
open model with its R matrix given by i?^ n_1 ^(w). Repeating the procedure j times, we 
can obtain the eigenvalue Fj(u^; {ti^}; {^2]}) with its transfer matrix constructed by 
R( n - j \u), 

ri(« W) ;{t>g-?};{t;g>}) 

+w U\u^\u^; {v^})A^(u^ {v%})r j+1 (u^; {t/g£>}) 
+w% j {u^% j {u^ ] {v^})A^(u^; {„$}) (69) 

with «W = « - 2jr,, = 4 J) - 277, {<>} = rf\ • • • {«} = R>}, 043} = 

{^mii} = {°}> m -i = N,m = m, 
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CO 



{fife 1 ?}) = {eg- 1 )}), 



(70) 



& V>; {<_?}) = n 



n 



a n -j{u 



a 



n-j 



1=1 Qn—j 



On-j[V t 



n-j 



m,-_i 



(U 



(?) 



^- 1) )e n .,(^)-^- 1) ; 



=1 a 



n-j 



(?) 



~0-i) 



)a„- 



(71) 



A^( u 0);{^)}) = g^-i-^-^+^On-j-i^')-^) 



A?(^ } ;{^}) = n 



e re _,(^)-#)e n _,(^) + ^) 



ic B _ i («0-)-^)c B _ i (ii(j) + ^) 

The coefficients iy's, w's vary with the corresponding transfer matrix. The Bethe anzatz 
equations are 

= l( °U)(^; { ^})Ar- 1 (^; {#})^ ^ (, = 1 '-'^ )(72) 

When j = n — 1, our problem become to obtain the eigenvalues of the nineteen- vertex 
models a£\ b{ 1] for 42, ^ and six- vertex models 4 2 \cf ) for 42-i,Cl 1) , respec- 
tively. For D^ 1 ) when j = n — 2, we will arrive at model whose eigenvalues can be 
given in terms of the product of the eigenvalues of two 4^ vertex models [8], [11], [21]. 
The six- vertex and nineteen- vertex models have been well solved [15], [17]. Thus, we can 
get the whole eigenvalues and the Bethe ansatz equations of transfer matrix for the 
42, 42-i, Bn\ CfP>D(p vertex model with open boundary conditions. 
From Eq.(63), we can obtain the energy of the Hamiltonian Eq.(12) 
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The more explicit expression for Eq.(73) is 

- ro ' (0) S co sh (, i -t)-l h (2,) -^ (0)1 - (74) 
2.5 The quantum-algebra-invariant cases 

Here we will consider the trivial K matrices for A^, A^-i, models (quantum- 

algebra-invariant case). The K matrices take the forms 

K- = l, K + = M(n). (75) 

We begin with A^, models. After a not very long derivation, we can rewrite Eq.(63) 
in detail, which is 

r(«) = r («) = ^ 0) (u (0) )4 V (0) )d°V (0) ; {«LV}M (Tno) («) 
+<H« (0) )^(« (0) )^ 0) (t* (0 M«fc 1 J})c (mo) («) 

n-2 

+ E ^(w (j V.(« 'V? +1) (« 0+1) )^ +1 V (j+1) )d 0) (« (0) ; {v { m ]l})B im - m ^\u) 

3=0 
n-2 

3=0 

+/i„_ 1 ( M ("- 1 ))^_ 1 ( M ("- 1 ))d° ) (« (0 M^L 1 1 ) })^ (m "- l) («), (76) 
where the identity A™ J (w (j) ; {v&]})$ +1 \uV +1 '>; {v<g}) = 1 is used, 

lij{u^) = f[ Uj{u^) = fl w {l \u {i) ), (77) 



fc=i e„(# - v\ + ') 

AA e n _,(^) + ^)) e „_ j( ^)-^)) 

X TT +Vi)On-j-l[ Vj_) 
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i\ e B _ J -_ 1 («0-+i) - ^ +1) )e n -,-i(^ +1) + tJp +1 >) ' 

C7 = 0,l,2,...,n-2) 

B (mn - l) (u)= TT 

where /o( M ) is just the matric elements f n (u) when n = 0. The more explicit expression 
of Bethe equations Eq.(72) are 



n 



On-j K + V k ) a n~j{Vi-Vi 



-, (i = l,...,m j; jVw-l)(82) 



n 



*=i &i(^- 1) +4- 2) )^(^- 1) -4 n - 2) ) 
m "" 1 ai(-i;{ n - 1) + ^ 1} )ai(-#" 1} " *f ei(if ~ 1} + ^^(if _1) " *f ~* , 



X II , („- 1 _f„_l) . („-J.; („-l), („_].) („_ |.) . („_J 

i^ihi-vl >+v} 'jhi-vl >-v\ ') f Q (vl > + v\ >)f {vl >-v\ ') 

-(n-l)^(n-l)^ 

= - _(„_!) („_!)' (n _ 1} , («= l,--,™n-i) (83) 
From Eqs.(39, 40, 41, 43) we can obtain 

= i. a »(««>) = 

sinh(ww; — 2?7) 



(?) 



e 2(2(n-j)-i)„ sin h( u (i)) cosh(7j^') - (2(n - j) + 3)r/) , 0/A 
l( m ) = 7777i777777i 777 773 ,,..u\ 77771 ^ , 77777 - I 84 ) 



2(n ^' )+lV ; sinh(uO') - 4(n - j)r/) cosh(wCj) - (2(ra - j) + l)ry) ' 

{j) {j) sinh(w^ - 2(2(ra - j) + l)r/) cosh(uk') - (2(ra - j) - l)rj) 
w i \ u ) 



sinh(-u^) — 2-q) cosh(-u^) — (2(n — j) + l)rf) 
ah(^) - 2(2(n + l)7y) 0) 
sinh(uC?) - 4(ti - j)rj) ' ™ 2 (^)+ l( 



= e- 2 -smh(^)- 2 (2(n-,) + 1)7,) ^ „ , )} = ^^.^ 
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and 



2e 2r > sinh(wp ) ) sinh(^ ; - 4(ra - j>) cosh(t^ - (2ra - 2j - l)r]) 



,(?) 



00 



sinh(f. 



sinh(u- n ^ — 2r/) 



^(u^) = 1, wW^'') 



sinh(up^ — 2rj) 



, j<n-2 

j = n-l 
(86) 



3^ sinh(w^) 



(J) 

2(n-j)+l 



sinh — 2?7) ' 
e 2( 2 (n-j)-i)r, s j n h( w O)) s inh(^') - (2(n - j) - 3)77) 
sinh(wtO — 4(n — j — 1)?]) sinh(w^) — (2(n — j) — l)rj) 



(87) 



(3) 



sinh(w^) - (2(n - j) + 1)?]) sinh(uk') - 2(2(n - j) - l)rj) 



w 



Cj)( u Cj)) = 



e 2,? sinh(it 



sinh(u^) — 2r/) sinh(it^) — (2(n — j) — 1)77) 

(?) 



2(2(n-j)-l) V ) U) 
sinh(u^) — 4(n — j — 



-2(2(n-j)-l 



>#8) 



2e 2,? sinh(v l °' ) ) smh(^ j) - 4(n - j - 1)?]) sinh(4 j) - (2n - 2j + 1)77) 



e 2ri smh.(v\ 



smh(i>^ — 2rj) 



, j<n-2 



sinh(vf n ^ — 2rf) 

for and respectively. The Eq.(63) for A^_ ± , and Z?W can be written as 

r(u) = v (u) = ^ V (0) )4 V {0) )d°V (0) ; R"_ 1 1 ) })^ (rno) («) 

+<) (t,<°>)u>(°> (ti<°>)^ (^ (0) ; })C (mo) («) 



j = n - 1 
(89) 



1 

n-2 



3=0 
n-2 



+ E ^(« 'V.(« (j VP +1 V ' +1) )^? +1) (« ' +1) )d 0) (« (0) ; {CA ) }) B(m, '' m ' +l) («) 

,C?+i) ^,0'+i)v-,0'+ 1 ) ^,(i+i)^(o)/ ,(0). r~(-i)- 



i=o 



and 



r(«) = r„(«) = 4V 0) )^ (0) )d°V 0) ; {t£l})A™ («) 
+<) («<°> )a>J» («<°> )d 0) (« (0) ; »C (mo) («) 



n— 3 



+ E ^(« 'Vi(« 'VP' +1 V (j+1) )^ +1) (« 0+1) )d V (0) ; {C_ 1 1 ) })^ (mj,m, ' +l) («) 

n— 3 



3=0 



u 



(i+i)v-,0'+i) 



u 
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+^-3(u (n - 3 Vn-3(u (n - 3 Vi n_1 V (n ~ 2) ^ 
+^S(« (B - S W3(« (B ^^ 



-2,m„-i 



respectively. Here Eqs. (80,81,82) hold until j = n — 3 for A^-i, and j = n — 4 
Di}\ The rest B's and B's are 



m»-2 



(«("- 2 )+i5ji n - 2) )ai(M (n - 2) -4 n " 2) ) 



x 



i e 2 (w(«- 2 ) + 4 n - 2) )e 2 («(«- 2 ) - v£~ 2) ) 



j g(m„_2,m„_i)^ = TT 1 I" ;.■ .« U" 



for 4^ 



gK-2,m„-i) ^ _ J~J 
fc=l 

x n 



fc =i e 2 («(»- 2 ) + 4 n - 2) )e 2 (fi(«- 2 ) - vt~ 2) ) 
X ^ ei^C"" 1 ) + v { r 1] )ei{u^) - v\ n - l) ) ' 

ai (u(n-2) + 5 (»-2)) ai ( i2 (n-2) _ 



fc =i e 2 (^" 2 ) + 4 n - 2) )e 2 (w(™- 2 ) - v ( r 2) ) 

a(_fi(n-i) + i)( n - 1 ))6(u(»-i) + v ( r 1] + 477) 



i=i fc(-« (n_1) + tf _1) )a(« (n_1) + *>i (n_1) + 



X 



fc =i e 2 («(™- 2 ) + v { r 2) )e 2 {u^) - ^ 2) ) 
m "" 1 a(u^ - v { r 1] )<u {n - l) + v ( r l) - 477) 



n 



for CM and 



i &(£(»-!) - j5, (n_1) )a(« (n - 1) + ui (n_1) + 477) ' 
3 a 2 (M( ra - 3 ) + v^ n -^)a 2 (u {n ' 3) - uf- 3 )) 



1 e 3 (M(™-3) + i#~ 3) )e 3 (7>-3) - ui n_3) ) 

m f T 2 a(-fi("- 2 ) + ^ n - 2) )a(-« (n " 2) - v ( r 2) ) 
x 11 



fe=i 

x n 

1=1 

m n -2 









" 2) )K- 


-«("- 


-2) 


~(n- 

-n 


- 2) ) 


a(- 




- 2) + *f 




_„(n- 


-2) 


~(n 
~ V l 






-7j(«' 






_7j(n- 


-2) 


~(n- 

-"1 




a(- 




- 2) + *T 


2 V(- 




■2) 


-<«- 
- n 


2) ) 



-Ji 6(-S(«- 2 ) + v { r 2) )b(-u("-V - vf 2h 
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m„_ 3 



C2 



(fc(n-3) + {)(«-3)) e2 (^(«-3) _ fl(n-3)) 



s 

rrin-2 

x n 



i e 3 (u(^ + #- 3) )e 3 (w(™- 3 ) - #" 3) ) 

a (£(n-2) + ^- 2 )) a (^(n-2) _ -(n- 2)) 



r,(™- 2 )w„^n-2l „~,("- 2 )^ 



gK-2,m»-i)( M ) = JJ 



" =1 6(m(™- 2 ) + 4" -">(m(™- 2 ) 

x \=i ^ (n_2) +v, (n_l) )6(«( n - 2 ) -v\ n - i] y 



m n - 2 



k-- 



\ & (fi(n-2) + 4«- 2 )) & (,i(n-2) _ £(- 2 )) 



x 'i b(-s(- 2 )-^ i) )6(-^- 2 )+{); n - i) )' 

(92) 

for where a(u) = sinh(|— 477), b(u) = sinh(|) for and a(u) = sinh(|— 2r/), 6(u) = 
sinh(|) for -D^ 1 - 1 . The other Bethe equations are 

m fr 3 q 2 fe (w - 2) + 4 n " 3 V 2 (^ ( "' 2) - vt 3) ) 

11 / ~(n-2) . ~(n-l)\ / ~(n-2) ~(n-l)\ 

m " T 2 a 2 (-^ ( "- 2) + i)(«- 2 ))a 2 (-{;("- 2) - 5(™- 2 )) e 2 (# " 2) + 6f - 2 ))e 2 (f;f ~ 2) - ^ 2 )) 
x 11 



% h(-vt 2) + &~ 2) M-vt 2) - vi n ~ 2) ) <vt 2) + v { r 2) ) ai {vt 2) - v { r 2) ) 

^- i) (i;f- 2 Vi(2^ 2) ))^^- 2 )(^- 2) )^r i) (#- 2) ) 

j l z — J-j ■ ' ' i m n-2) \"'->) 



fc=i 

= l,...,77l n _i) (94) 



X 



^ 2 
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for 4ti, 



~(n-2) ~(n-3) 



)a 2 (6. 



~(n-2) ~(n-3) 



~(n-2) _ ~(n-3)x 
k , 



m n -i 1/ ~(n— 2) ~(n— 1)\ /~(n— 2) ~(n— 1) . , \ 

x -q K- V l + v i >K >+4r]) 



1=1 



a(-v\ 



~(n-2) ~(n-l)\ 7 /~(n-2) ~(n-l) 



7J 



+ 4?y) 



"ln-2 



X 



a 2 



(-# " 2) + ^r 2) )a 2 (-vt~ 2) - i { r 2) ) e 2 (^ n ~ 2) + ui n - 2) )e 2 (if " 2) - i^" 2 )) 



-(n-2) 



-(n-2) . ~(n-2) 



w (-l)(^-2) )ai(2 ^-2) )) ^(n-2) (u (n-2) ) -(n-l) ( ~(n-2) ) 



\ /~(n— 2) ~m— zi\ 

)ai(u> '-vi ; ) 
1,- ■ • ,rn n - 2 ) (95) 



/?n-lK (n - 2) ) 



_(n-2) /(i ,(n-2)> 



'(«T") 



ai({;, 



m„_ 2 
II I 



~(n-l) . ~(n-2) 



_(n-l) ~(n-2) 



- VI 



+ vt 2) 



)e 1 (v. 



~(n-l) ~(n-2) 



x n a( ' 



~(n-l) ~(n-l) 



')&(-«. 



~(n-l) 



-47]) 



/~(n-l) 

aw; 



)a(5 



~(n-l) . ~(n-l) 



-(n-1)/ (n-1), 
^2 K ') 



(n-1)/ (n-1) 



')A»(u 



(n-1), 



(i = l, 



-4?7) 

A-i) 



(96) 



for C« and 

m n _4 

n 



fc=i ^ 
x n 

Z=l 

m„_i 

x n 

t=i 



/~(n-3) 




- 3 ) _ 


~(n- 


- 4) ) 




:(^" 3) 




-3) _ 


~ (n- 


■ 4) ) 




fe(-# 


-3) . ~(n-2) w 


~(n- 


-3) _ 


it 




/ ~(n 


-3) . ~(n-2)x / 

>+v\ >(- 


~(n- 


-3) _ 


~ (n- 

vi 


- 2) ) 


b(-t 


-3) . ~(n-l) w 

+ vl >)b(- 


~(n- 
u i 


-3) _ 


~(n- 
Vt 


"") 


a(-v\ 


-3) . ~(ra-l)\ , 

+ vt )°(- 


~(n- 


-3) _ 


~(n- 

v t 





x "g 3 a 3 (-^ (n ' 3) + ^ n ' 3) )a 3 (-^ ( "' 3) - ^"- 3) ) e 3 (#" 3) + ^"~ 3) )e 3 (if " 3) - v^) 



1% h(-vt 3) + vi n - 3) )b 3 (-v. 



(n-3) 



-(n-3)x 



Vs ')a 2 (v ( [ l 3) + t4™ 3) )a 2 (^, 



(n-3) 



^ 3) 



) 



(i 



■ , m n - 3 ) 



(n-3)/ (n-3)\ 

K ) 



(97) 



mn-3 /~(n-2) ~(n-3)\ /~(n-2) ~(n-3) 



n 

fc=i 



a v 



+ v 



a f ; 



_(n-2) 



+ 4^ 3) 



)b(v. 



(n-2) ~(n-3)N 
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X 



n 



(n-2) ~(n-2^ 



& a(^-^- 2 ))a(^- 2 ) + ^- 2) ) 



^i(4 n - 2) )^-i(^r zj ) 



(98) 



n 



X 



i 6(uJ n_1) + 4 n_3) )6(f;f _1) - 4 n_3) ) 

m„_i ^ ~(n-l) , ~(n-l)^^ ~(n-l) _ ~(n-l)> 
II 



a(-v> '+v\ ')a{-v. 



Q?(4 n - 1) )/3n-l(^ n " l; ) 



(n-1) 



-. (i = 1,. • - ,m n _i) 



for £>W. 

The coefficients are 



(99) 



(li^) = 1, U®(u®) 



e 2v sinh(it^ 



sinh(wO') - 2?]) ' 
e 2(2(n-j))» sin h( M (j)) C0S h( M (3) - 2((n - j) - 1)77) 
sinh(M^) - 2(2(n - j) - 1)77) cosh(u^) - 2(n - .7)77) 



(100) 



w <j\ u W) = S -^ uU) ~ 2(2(n " j))7]) cosh ^ 0) - 2((n - j) + l)r/) 



w (j)( u Cj)) 



sinh(wO') - 2?y) cosh(itO') - 2(n - j)rj) 

, ^L)^ ) = e^ 2 ^'", (iOl) 



e" 2?? sinh(^) - 2(2(n - j))??) 
sinh( M 0) -2(2(n- j) - 1)77) ' " 2( 1 



2e 2tl sinh^) sinh(^ ; - 2(2(n - j) - l)rj) cosh(u^ ; - 2((ra - j) + 1)77) 



.0") 



(j) 



e 4 ^ sinh(2u> 



(n-l)x 



sinh(2vf n 1 - ) — An) 



for A 



(2) 

2n-l) 



sinh(i^ — 2 77) 



j < n — 2 
j = n-1 



(102) 



= 1, u {j \u^) 



e 2r > sinh(w^) 



sinh(wO') - 277) ' 
e 2(2(n-j)) 7? sin h( w 0)) sinh(M^) - 2((n - j) + 2)77) 
mnh(uO') - 2(2(n - j) + 1)77) sinh^C?) - 2((n - j) + 1)77) 



(103) 
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mnh(«w) - 2(2(n - j) + 2)77) sinh(?/W - 2(n - j» 

sinh(wO') - 27/) siiih(uO') - 2((n - j) + 1)77) 
e-^sinh(^)-2(2(n-j) + 2)r ? ) 0) W) 



-2(2(n-j)) V 



for C« and 



sinh(uO') - 2(2(n - j) + l)r/) 
2e 2,? sinh(4 i} ) sinh(^ j) - 2(2(n - j) + 1)77) sinh(v l (i) - 2(n - j)rj) 



(104) 



e 4l) sinh(f|" 1 - ) ) 
sinh(v| n-1 ^ — 4r/) 



sinh(up' ) - 2rj) 



, j<n-2 

j = n-l 
(105) 



with 



e 2r > mnh(u^) 



0> 



^ )(M0)) "smh(^)-2r/)' 

e 2(2(n-j)-2)T, s i n h( W 0')) sillh(u^) - 2((n - j) - 2)7/) 

sinh(utf) - 2(2(n - j) - 3)77) sinh(uk') - 2((n - j) - l)r]) 

(j<n- 3) (106) 



Cj)( u (j-)) 
w 0)( u 0)) 



sinh(w^) - 2(2(n - j) - 2)77) sinh(u^) - 2(n - j» 



sinh(w^) — 2rj) sinh(uC?) — 2((n — j) — 3)r/) 
" 2 " sinh(w^ - 2(2(n - j) - 2)77) _ {j) 



sinh(uO') -2(2(n-j) -3)77) 



W 2(n-,)^ 



00) 



-2(2(n-j)-2)j, 



(j<n- 3) (107) 



w (»-2)( u (n-2)^(n-2) (u( ^ 

2 ))^ n - 1) ( M (n ~ 2) ) = 2^ 2 1 ( M ("- 2 ))^ 1 2 ( M ("- 2 ))^( M ("- 2 ))fi 2 ( M ("- 2 )), 
2^ 1 1 (^- 2 ))^ 2 ( M ("- 2 ))fi}( W ("- 2 ))fi 2 ( M ("" 2 )), 



(n-l)/ ( T 



("- 2 )^,(«- 2 ))-,("- 2 )^,("- 2 ) 



4 n - 1) ( W («- 2 ))^ n - 1) ( W ^ 2 )) = 2^ 2 1 ( W ("- 2 ))^ 2 2 ( M ("- 2 ))^( M ("- 2 ))fi 2 ( M ("™ 2 )) 



n}(u) = nl(u) = 1, 

sinh(-u — 4t/) 



fi 2 ( W )=ft 2 ( M ) 



e 2r ? sinh(w) 
sinh(-u — 277) 



sinh(tt — 277) ' 



Wftu) = («) = e~ 2 ". 



(108) 



2e 2ri sinh(?j l ' ) ) sinh(v^ - 2(2(n - j) - 3)t/) sinh(^ j - 2(ri - ,7)77) 



,00 



00 



sinh(f l - 



sinh(w^ — 2t/) 



sinh(t> P — 27/) 

j = n — 2, n — 1 



j < n - 3 



(109) 
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for Up to now, we have gotten the whole eigenvalues and the Bethe equations of 

transfer matrix for the A^, A^l-i, D^> vertex model with trivial K matrix. 

We have checked that our results can agree with that obtained by the analytical Bethe 
ansatz[20]. For the non-trivial K matrices, we still can repeat the similar procedure as we 
have done above. In the following section we will present some results for the non-trivial 
cases. 



2.6 Results for B^\C^\D^ models with diagonal K-matrices 
(one free parameter case) 

In this section we will consider the B^ , , models with K matrices having one free 
parameter. We do not discuss other non-trivial cases for the sake of simplicity. We can 
find the expressions of eigenvalues Eqs. (76,90,91) and Bethe equations Eqs. (82,83,95-99) 
still hold for , , models. So we only present the coefficients for these models. 
For B^\ from Eq.(17), we take 

K®(u, n, C+) = K { *\-u - n, C+)*M(n). (110) 
Then we can obtain the coefficients 

0}V°>) = e-«sinh(C- + |) sinh(C_ + \ - (2n - %) 

-(0), (cm sinh(^) smh(C- - f + 2ry) sinh(C- + f - (2n - 3) V ) 

U [U } ~ sinh( M -2r/) ' 1 j 

tf£U« (0) ) = e u sinh(C_ - |) sinh(C- ~\-{2n- Z)rj) 

2e u sinh((2n - l)rj) sinh(2r]) sinh(C_ + f ) sinh(C_ + f - (2n - 3)77) 

sinh(w — (2n — l)f]) sinh(« — 4(n — l)r/) 
e M sinh((2n - 1)77) cosh((2n - 3)77) sinh(C_ - f ) sinh(C_ + f - (2ra - 3)r?) 

sinh(w — 4(n — 1)77) 

U7 W( tl (0)) = flUsiT,W "- r9n+1 ^ 



( o) (0) _ e u sinh(w - (2tz + 1)77) sinh(7j - 2(2n - 1)77) sinh(C+ - f ) sinh(C+ - f - (2ra - 3)77) 



sinh(-u — 277) sinh(-u — (2n — 1)77) 
(0) (0) = sinh(^ - 2(2n - 1)77) sinh(C + - f ) sinh(C + + | - (2n - 1)77) 

1 J sinh( M - 4(n - 1)77) ' 1 ' 

«>2U« (0) ) = e"" sinh(C + + £ - (2n - 1)77) sinh(C + + £ - 4(tz - 1)77), (113) 



2 v , „ 2 



(o) 2e Vl sinh(t>j) sinh(t> j — 4(n — 1)77) sinh^ — (2n + 1)77) sinh(C+ — y — (2ra — 3)77) 



sinh(t>j - 277) sinh(C+ + y - (2n - 1)77) 

(114) 

While for 1 < j < n - 1, Eqs. (87-89) still hold. 
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Similarly, for Cg\ by Eqs.(18) and (110), we have 

^V>) = « " smh( ( _« + — ), „(.)(„<.)) = 



(115) 



e 2(2(n-j)-i)TH- V sinh(u^)) sinh(u^) - 2((n - j) + 2)r]) 
sinh(wO') - 2(2(n - j) + 1)77) sinh(ifCJ') - 2((n - j) + 1)77) 

x sinh(C 0) - — + 2((n - j) + 1)77). j < n - 1 



e^ 1 sinh(Ci j) - ^-) sinh(^) - 2(2(n - j) + 2)77) sinh(u^) - 2(n - j» 



sinh^-?) — 277) sinh(ii(-7') — 2((n — j) + 1)77) 



w 



(?)(„(*)) 



e - " sinh(u w - 2(2(ra - j) + 2)77) 
inh(uk') -2(2(tz -j) + 1)77) ' 



w. 



2(n-j) 



sin 

e 2 



j <n-l 



-2((n-j)-l>, ;„^ A (j) M 



0) 



'sinh(C; J; + 



2((n-j) + l)77), 



(116) 



2e 2^ ) +r ) S j n h( v ^) s inh(^ Uj - 2(2(n - j) + 1)77) 



,0) 



smh(v^ — 277) 



x sinh(^ j - 2(ra - j)rj) sinh(C+ ) - ! %-), j<n-2 (117) 

(n 



„0) 



e * smh(v- ; ) smh(C 



sinh(4 n - 1} - 477) sinh(CT" i; + 



(n-l) « 



(n-l) 



477) 



j = n-l 



where = ( T ± ji] 

For D^, there are seven 1-parameter K matrices. For the sake of simplicity, we 
only consider one K matrix. We choose Eq.(25) and the corresponding K+\u,n,£+) by 
Eq.(llO). Then we get 

uPiu®) = e- uU) sinh(Ci j) + ^) sinh(C a) - ^ + (2(w - j) - 4)77), 



a; 



Cj)( u Cj)) 



sin 



h(«(J)) 



sinh(wO') - 277) ' 

e uU) sinh(TjW) S mh(u ( ^ - 2((n - j) - 2)77) 



(118) 



sinh(wO') - 2(2(n - j) - 3)r/) sinh^tO - 2((ra - j) - 1)77) 



x sinh(Ci j) + 2t?) sinh(C W + (2(n - j) - 1)77) (j < n - 3) 

2 2 



{j) {j) _ e uU) sinh(M^') - 2(n - j)rj) sinh(7jW - 4((n - j) - 1)77) 



sinh(-u( J ') — 277) sinh(7j(-?) — 2((n — j) — 1)77) 
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x sinh(£^ + 



u 



sinh(w( J ) — 4((n — j) 



+ 

i)n) 



u 



(i) 



(j) 

2 , sinh(Cy ; - 2 



+ (2(n-j)-4)rj), 



sinh(uO') - 2(2( 



J) " 3)17) ' 

CO 

0')) = e -« (i) S inh(ci j) + 2rj) sinh(C 

2 



n 



(119) 



(?) 



0) 



+ 2((n-j)-l)rj), (j < n - 3) 



^( M (-2)) = e - 2 V 1 sinh(Ci n - 2) + 



7/, 



(n-2) 



,(n-2) 



^(t^^) 



e" 2 sinh(w^ n 2 ))sinh((_ 



„(n-2) 



+ 27/) 



sinh(«( n - 2 ) - 2r/) 



W^u^" 2 )) = 



e 2 



sinh(-u(™ 2 ) — 4r/) sinh((. 



(n-2) „(n-2) 



sinh(w( n - 2 ) - 2rj) 



^(u^) = e"^ sinh(Cr 2) + ^ - 277), 



fi 2 ( w (- 2 )) = e-T sin h(ci n - 2) 



2 

M (n-2) 



), 



u (""2) 



fi 2 ( W ("- 2 )) 



sinh( W ("- 2 )) sinh(cl" _2) + - 2r/) 



sinh(-u( n 2 ) — 2r/) 



e 2 



sinh(-u(™ 2 ) — 4r/) sinh(QT ' + 



(n-2) u (n-2) . 



sinh(i(( n 2 ) — 277) 



(n-2) 



iy|( M ("- 2 )) = e" 1 ^ sinh(C 



(n-2) 



7J 



(n-2) 



+ 277), 



(120) 



U) 



2e 2v "' sinh(^ (i) ) sinh(^ UJ - 2(2(n - j) - 3)r]) sinh(^ UJ - 2(n - j)rj) 



.(3) 



(j) 



sinh(i;^ — 2rj) 



x sinh(C? + ^) sinh(d j) - ^ + (2(n - j) - 4)77), j<n-3 
smh(v- ; ) smh(C+ ^ — ) 



smh(4 n " 2) - 277) sinh(d n ' 2) + - 277) 

smh(7j l > ; ) sinh(Q_ + - L "2— 

sinh(^ n_1) - 277) sinh(C| n ~ 2) - ^) + 277) 



, j = n - 2 



-, j = n - 1 



(121) 

where = C=f + J 7 ?- Now we present some ABA solutions to open vertex 
models with non-trivial diagonal K matrices(one free parameter case). Repeating the 
same program demonstrated in sections(2.2-2.4), we can obtain the ABA solutions for 
other non-trivial cases. 
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3 Conclusions 



In the framework of algebraic Bethe ansatz, we solve the A^, A%^-i, ver- 
tex model with trivial K matrix and find that our results agree with that obtained by 
analytical Bethe ansatz method [20]. Some results are presented for the C£\ 
with their diagonal K matrices having one free parameter. Using these results we can 
study the boundary effects etc. 

Here some open vertex models, such as are not considered due to their com- 

plicated R matrices. For these models, we need to generalize the ABA further. Addi- 
tionally, for the ABA has been generalized to the spin chain with non-diagonal reflecting 
matrices[27]-[30], it is an interesting problem on how to apply the method to other models 
with higher rank algebras. 
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A R matric elements 

The R matric elements are 

a n (u) 

b n (u) 
c„(u,i) 
c n (u,i) 
dju,i,j) 
d n (u,i,j) 
e n (u) 

fn(u) 
9n(u) 

9n( U ) 

forflW.tfW.DW and 

a n (u) 
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= 2 sinh(- - 2rj) sinh(- - K rj), 
= 2 smn(— J smh(— — kt]), 

Z Z 

= d n (u, i, i) — 2e~i sinh(2r]) sinh(— — kij), 

Z 

= d n (u, i, i) — 2e2 sinh(2r]) sinh(— — kt]), 

z 

= 2e l e j e-^ K+2 ^ ))r > sinh(2r/) sinh(-), 

z 

= 2e l e j e^ +(2(l - ]) - K)v smh(2rj) sinh(-) , 

z\ 

= 2sinh(-)sinh(--(/c-2)77), 

= b n {u) + 2 sinh(2r]) sinh(Kry), 

= — 2e~2 sinh(2?7) sinh(^ — nrj), 

z 

u Vb 

= — 2e2 sinh(2?7) sinh( kt\) (A.1) 

2 

. . ,u n . , ,u . 
= 2smh( 2r]) cosh( nrj), 



b n (u) = 2sinh(|) cosh(^ - icq), 

c n (u,i) = d n (u, i, i) — 2e~z sinh(2r^) cosh(— — nrf), 

c n (u,i) = d n (u, i, i) — 2e2 sinh(2^) cosh(— — nrj), 

d n (u,i,j) = -2e J e j e-t + ( K+ ^-^sinh(2r])sinh(^), 

d n (u,i,j) = 2e i e J e5 + ( 2 ^-^sinh(27 7 )smh(|), 

u u 

e n (u) = 2sinh(-)cosh(- -(k-2)t]), 
fn( u ) — b n (u) — 2 smh(2r]) cosh^KT]) , 
9n{u) = -2e~^ sinh(27]) cosh(- - nrj), 

g n {u) = -2e% sinh(2?7) cosh(- - nrj) (A.2) 

for 4«, 4S-i- Where a = 1 for 42, £«, models, Ci = 1 if 1 < i < n and a = -1 
if n + 1 < i < 2n for 42-1, models. 

B Derivations of diagonal operators acting on eigen- 
vector 

Acting the diagonal operators x(u) = A(u), A aa (u), A 2 (u) on the m-particle state and 
having carried out a very involved analysis as that in Ref. [15], we can obtain the following 
expression 

x{u)\T m {vi, - ■ ■ ,v m )) = \^ x (u,{v m })) 

m 

+ E h l ( u , v h <0|*m-l( u > ^ {^m})dd) 
i=l 
m 

+ E k 2( U ' V ^ ^l^m-lK V i'> { V m})dd} 
i=l 
m 

+ J2K{u,Vi,a x )\^_i{u,Vi; {v m }) 



i=l 
m 

+ E ^(^ V i> «x)l*S-l(«, ^ {"m}) 
i=l 

m—1 m 

+ E E ^ e2 ^M 1 ( M '^' V i)l^™-2(«,^ ) ^;{Vm})d 1 e 2 ) 
i=l j=i+l 

m—1 m 

+ E E ^K^^.OSdJ^UK^,^;^™})^/ 

i=l j=i+l 
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m— 1 m 

i=l j=i+l 
m 1 m 

+ E E ^(W, «i)/di l*m-2(«, Vi, Vf, {V m }) e d d lf ), 
i=l j=i+l 

where when x = A, A aa , A 2 , 

h*(u, Vi, d) = a\(u, Vi), R?(u, vi)% a\{u, Vi), 
h%(u, Vi, d) = al(u, v^, R^(u, Vi) d a a d , a\{u, v t ), 
h%(u, v h a x ) = 0, Rf(u, Vi, a), a\(u, v h d), 
hl(u, Vi, a x ) = 0, Rf(u, Vi, a), a\(u, v u d), 
a x = 0, a, d, 

' Wi(m)A^(u; v u ■ ■ ■ , v m )\T m (vi, v m )), 

|**(«,R»})> = ^••■^(^.••^«m)[T^(«;{t; m })J:^]aai^"- 6 "|0>, 

k u q (u)A™(u; v ir --, v m )\T m (vi, v m )), 
respectively, and we denote 

l*m-l(«> ^ K»})/ dl > = B/^) • • • , Vi, ■ ■ ■ ,V m ) 

xsf-tivi, {^})Ar 1 K; {*i})viMF*"*"\0), 
Kl) /«*) = £/H$m-T m K • • • , ^, • • • , U m ) 



|*m-2(«, ^ {Vm})d ie2 ) = F(«)^|™(t;i, • • • , U i} • • • , Vj, • • • , 



x [f m -\ Vl ; {Vi, v^ZZZheSZ'TJvr, ^i, v^Sftivi; {Vi}) 

xAr 2 K;{^,^})AK)F 6 -- 6 -|0), 
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\^ ( J- 2 (u,Vi,Vj; {v m }) fdl ) = F(u)§Z-2 m {vi, ■■■,v i ,---,v j ,---,v m ) 

xAr 2 (^;{^,^})A(^)F bl - 6 -|0), (B.9) 

|*S- 2 («, "i, ^ K,})*/) = • • • , ^, • • • , Vj, ■ ■ ■ , V m ) 

x \f m -\ Vl ; {v t , v 3 })zZzUf m -\v 3 - { Vi , r,}t::Z j r , 

xsZ'Zi^ i*i,*j})s£::£(v i ; w^-Mo). (b.io) 

The explicit expressions of Hf dl (u, v i} Vj), I — 1,2,3, A are listed as below 

H l,di( U > V ii V j) = a l( U > V " dl){c\(Vi, Vj) + c)(Vi, Vj)) 

+a\{u, Vi)(cl(vi, Vj, di) + cl(vi, Vj, d ± )) + b\(u, v^g^Vi, Vj, d x ) 

+a\(u,v i )a\(u,Vj)g 1 (u,v i ,d)r(v i - u) d d di , (B.ll) 

H^ di (u, v u Vj){ e 2dl = a\{u, Vi, di)(cj(ui, Vj) + c\{vi, Vj))S dlf 
+al(u,v i )(R^(v i ,Vj) f c ^ dl + R%(vi,Vj) f c * dl ) 

+b\(u, Vi)gi(vi, Vj, di)5 dl c 2 + a\{u, Vi)a\{u, Vj)g 2 (u, v, t , f)r(v t - u) f C2dl , (B.12) 

H 3,di( u ^ v ^ v j) = aKu^^d^cKv^Vj) +al(u,v i )c*(v i ,Vj,d 1 ) 

+b\(u,v i )g 2 {vi,Vj,d 1 ) + a\(u,v i )a\(u,Vj)g 1 (u,Vi,d)r(v i - u)f idi , (B.13) 

H ldA u i v ^ v i)fdi = a l( u ' v ^ di)c\ (vi, Vj)8 dld + a\(u, Vi)R%(vi, Vj)f dl 

+bl(u, Vi)g 2 (vi, Vj, di)8 dlf + a\(u, Vi)a\(u, Vj)g 2 (u, v h d)r(vi - u)f di , (B.14) 

H i,di( u i v i> v j) = Riiu^'dtai^i^j) + c \{ v h v i)) + b 2 2 (u,v i )g 1 (vi,Vj,d 1 ) 

+Rq(u, Vi)(c\(vi, Vj, d x ) + cl(vi, Vj, d ± )) + f(u + Vi) d d f(u - Vi) d dl a R^(u, Vj, Ji)e}(^, u, d) 
+f(u + v t y d e g f(u - ViY^Mu, h ) R i(^ Vj)f if f( Vi - u)^, (B.15) 

H*%(u, Vi, Vj) f c e 2dl = Ri(u, Vi)%M(vi, Vj) + cl( Vi , Vj))8 dlf 

+Rq(u, Vi){R^(vi, Vj){ e 2dl + R%(vi, Vj) f C2dl ) + R((u, ViY^g^Vi, Vj, d^S^ 
+r(u + v t )T g f(u - v^Mu, v h f)R?(u, v^r^ - u) f b e c 

+r(u + Vi)? g f(u - Vi)ZR£(u, Vj, c 2 )Rl e (v t , «)£, (B.16) 

Hf%(u, Vi, Vj) = Rf(u, Vi) a d d a c\(vi, Vj) + R£(u, Vi)c 2 7 (vi, Vj , Ji) 

+b 2 2 (u, Vi)g 2 (vi, Vj, d^ + r(u + Vi) a d d f(u - Vi) a d dl a Rf(u, Vj, di)e{(vi, u, d) 

+f(u + Vi )T g r(u - ViYlM^ v h h)R^(u, Vj)f if f(vi - (B.17) 
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+i?f («, vJStoivi, Vj , d^j + f(u + vX c g r(u - v t y d ( a Ri(u, Vj , f)R be (v t , u)f c 

+r(u + vX 9 r{u ~ v^tM^ d)R*{u, v 3 )f h f(v t - «)* (B.18) 

^K^^i) = a^(M,fi, Ji)(c5(^, ^) + c4(ui,Vj)) + 

+a 3 (u, Vi)(cl(vi, vj, di) + cl(vi, vj, di)) + a^w, Vi)a\(u, Vj, di)e\(v h u, d x ) 
+a 3 1 (u,Vi)al(u,v j )g 1 (u,v i ,d)r(v i - u)f idi , (B.19) 

H 2,k K V i> V j) f C2di = a e( M ' V i> di){c\{v U Vj) + cl(Vi, Vj))8 dlf 

+a 3 (u,v i )(R^(v i ,v j ) f c ^ dl + R^v^v^J 

■,Vj,c 2 )R b 3 e (v i ,u Jc2 , 

H^dM^i^j) = (4,{ u -> v i-,di)cl{v h Vj) +a%(u,v i )($(v i ,v j ,d 1 ) 
+b\(u, Vi)g 2 (vi, Vj, di) + a\(u, v^)a\(u, Vj, d 1 )e\{v i , u, d L ) 

+a\{u,v i )al(u,Vj)g 1 {u,v i ,d)r(v l - u)f di , (B.21) 

H* 2 di (u, v h Vj)f di = a 3 (u, v h d^cl^Vi, Vj)5 dld + af (u, v^R^Vi, Vj)f di 
+bl(u, Vi)g 2 (vi, Vj, d l )5 d j + a\(u, Vi)a 3 5 (u, Vj, f)R b 3 (vi, u)f dl 

+al(u,v i )al(u,Vj)g 2 (u,v i ,d)r(v i - u)f dl . (B.22) 
In Eqs.(A.ll)-(A.22), the repeated indices a and d\ do not sum. We can check that 

H^{u,v u vi)% a Hl di (u,v t , Vj ) d d l d d \ 



+b 3 (u, Vi)gi{vi, Vj, d 1 )5 dlS2 + a\(u, Vi)a 3 A (u, Vj, c 2 )R b 3 e (vi, u){ e 
+a\(u,Vi)a 3 2 (u,Vj)g 2 (u,Vi, f)r(vi - u) f c * di , (B.20) 



(B.23) 



(B.24) 



In Eqs. (B.23, B.24), all the repeated indices do not sum. We conclude that Eq.(B.l) 
can be verified directly by using mathematical induction, although it is a rather hard 
work. Similar to assumption of algebraic Bethe ansatz, we might assume that "quasi" 
m-particle states such as 5$ m _i|0), £$ m _i|0), BB$ m _ 2 \0), BE$ m _ 2 \0), EB® m _ 2 \0), 
F$ m _ 2 |0), FB& m _ 3 \0) etc are linearly independent. Here all the indices are omitted and 
all the spectrum parameters in the "quasi" n-particle state keep the order {v^ , v j 2 , • • • , v i k } 
with zi < i 2 < ■ ■ ■ < i k . For example, Bi(vi)Bi(v 2 )§ b £Z b 2 m (v 3 , ■ ■ ■ ,v m )F ub3 - bm \0) and 
B 1 (ui) J B 2 (w2)3 , m-2 m ( t '3' • " ' v m )F 12b3 - bm \0) are thought to be linearly independent. Then, 
by using this assumption, the property of Eq.(58) and some necessary relations, we can 
prove the conclusions Eq.(B.l) as done in Ref.[15]. 
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In order to obtain the eigenvalue and the corresponding Bethe equations, we need to 
change Eq.(B.l) by the following four useful relations (the proofs are omitted here) , 

sS 1 ^r(«i;W)ri(fi i ;{t; m })g::s: = 

(pLm-'T^iv^T^vf, m^^sti^ m, (B.25) 



x[T m ->,; {v^v,})^;^^:^ {v^stlZt:^ {*}), (B.26) 

sZ-Z^v {*i,tj})sz:.Z(^ «M«,-; {v m })tt = 

AT m ~ 2 {v,- {WMXZ^^^v* {v^sttiv^ {*}), (B.27) 

*z:;;:(r.r R v^s^:::^ MM*; {*J)r c £M%; {v m })tt 

= (p n -i(vi - v 3 )p n _M)- 1 T^\vdT^\v 3 )R^\v % + v^R^ivj + 
xStjriv,; {v^stlZt^ {*}), (B.28) 



where 



ci---c m J 

= m ■ ■ ■ L ^ v m tz-^ k-{u) hm 

*L-\-u,v m ))2Z m L-\-u^^^^ (B.29) 
T^\ Vl ) = kj(v, i )R^- 1 \2v l - 4r/)S. (B.30) 



Let 



A™(«; {v m }) = ]J Pn-i(u - Vi)p(u, v^, (B.31) 

i=l 

Pn-i(u) = a n _i(w)a n _i(-M), p{u,v) = — - — ■ — ^— -, (B.32) 

a n (u + v)e n {u - v) 



31 



v t ; {v m }) fd ) = ^S^ Ar 1 ^; WW^^to, • • • , • • • , vm) 

xSfZ m (vi; {Vi})n(vi, {v m })tt F b ^\0), (B.33) 
i 

l*£Li(«, « 4 ; K}U> = ^r^y ^" 1 ^; {^})^(«)Ci Cm (^i, • • • , • • • , v m ) 

\¥*l 2 (u, Vi , Vj] {v m }) dl ) = F{u)^^(v u ■■■,v i ,---,v j ,---,v m ) 

xA?- 1 ^; {^})Ar 1 K-; {v^UivM^F^lO), (B.35) 
\^_ 2 (u,Vi,Vj; {v m }) dl ) = F(u)$ e ^2 m {vi, ■ ■ ■ ,Vi, ■ ■ ■ ,vj, ■ ■ ■ , v m ) 

xst 3 Z m (^ {t^sttivu mri^ {v m }) c b \::z 

xAr 1 ^! {v t })u{v t )ptMK~\^ {^•})o; 1 (^)F 6l - 6 -|0) ) (B.36) 
\¥^_ 2 (u, v h Vj ; {v m }) dl ) = F( M )$^: 2 em (v 1 , ■■■,v i ,---,v j ,---,v m ) 

x^Tte; {v h v 3 })s d c -::t(v t ; M)^-; {v m })t\;; c c 

xAr 1 ^; {v i })u(v i )pi_ x $)\r 1 i?i\ {vi})^{vi)F^"^\% (B.37) 

|*m-2(«, u i> v j\ {Vm})^) = F(u)^Z 2 m (vi , • • • , U i} • • • , Vj, ■ ■ ■ , V m ) 

xs^z-ivr, { Vt , M)^; {v m }) c a \:: c z 

xnivr, {vMlZlZK- 1 ^, {^})A-- 1 K-; {v^coivM^Pn-mF^lO). (B.38) 

Using relation Eq.(B.25), we can easily change the \^m-i{ u i v i\ {"m}) jd) and vf, {v 
in Eq.(B.l) into |^J?-iK v u {v m })fd) and vf, {v m }) ab ), respectively. 
Let e 2 = Ci, from Eq.(B.26), we can get 

R^'Hh - v^t™- 2 ^-, {vi^y^u^;:^ {vi,^}^:^-, { Vl }) 

Pn-l(Vi - ^)ptl(O) R^K-Vi - Vj)T 
Pn-l(Vi+Vj) T^(Vi) 

x*t:Z:(<\r ^i, vA)Stt m (^ mMw, {v m })tt- (B-39) 

One should note that the repeated indices C\ and a,\ in Eqs.(B.39, B.42, B.45) do not sum. 
Using Eq.(B.39) and Eq.(B.23), we have 

H 2,di ( M ' V ^ V i) f cld, l*S-2(«» ^ { V m})d lC2 ) = 

u,Vi,Vj;{v m }) dl ), (B.40) 
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where 



Let e2 = ci, from Eq.(B.27), we obtain 

x^fcrrfe; {^,^})^. 2 i m (^; {vi})}ri(v j; {v m })ti- ( B - 42 ) 

Then, by Eq.(B.42), we have 

H tdMi V ^ V j)\^Z-2(u, Vi, Vj\ {V m }) dldl ) = 

Hdi( u > v i> v i)\*m-2(u,Vi,Vj; {v m }) dl ), (B.43) 

where 



HZ, (u, vA = ^ \ / . 1 ,J'T\ ■ B.44 

3 '" lV 3> pivj^JaKv^v^p^Vj - vjp^ivj + v i )T&(v j ) y J 

Let e 2 = oi, from Eq.(B.28), we achieve 

Pn—l\ u i ~r fjj 

x .... 1 t ' ] ' h ' St ■■;;■:(,■,: {vi,v 3 })S^ g y(v f , {*<})] 

\d\---dr, 



TW(vi)Tfa)(vj 



XT1(W, {v m })ZT m Tl(vf, {vMlZt- (B-45) 

Using Eq.(B.45) and considering Eq.(B.24), we have 

H% dl (u, Vi , Vj )f dl \*®_ 2 (u, v h v 3 ; {v m }) e d d lf ) = 

u,Vi,Vj;{v m }) dl ), (B.46) 

where 



Hl di (u,Vi, Vj ) 



piv^V^piVjjV^Pn-^Vj -Vi)p n -l(Vi + Vj) 



X 



Hf d {u, Vi, Vj) dd 



fl(n-l)(_£ _ fi.yd 



didi 



(B.47) 
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We can easily get 

S d ie2 Hl dl (u, Vi , Vj)\¥®_ 2 (u, vt, {v m }) die2 ) = 

(u,Vi,Vj;{v m }) dl ), (B.48) 



where 

H ld {u,v i ,v j ) = ^ rr ! — - T? r. B.49 

After making the notation 

it r r iu \ \^x{u,{v m })),x = A,A 2 

\^x(u, {v m })} = < J J _ (B.50) 

\ u(u)A?(u; {v m })^- dm (vi, ■ ■ -,v m )[T m (u- {v m })^\ aa F^^\Q),x = A aa 

we arrive at the final result Eq.(59). 
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